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Coulomb blockade thermometry using a 
two-dimensional array of tunnel junctions 

Tobias Bergsten, Tord Claeson and Per Delsing 

Department of Microelectronics and Nanoscience, Chalmers University of Technology 
and Goteborg University, SE-412 96 Goteborg, Suede 



We have measured current- voltage characteristics of two-dimensional arrays of small tunnel junctions 
at temperatures from 1.5 K to 4.2 K. This corresponds to thermal energies larger than the charging 
energy. We show that 2D-arrays can be used as primary thermometers in the same way as ID-arrays, 
and even have some advantages over ID-arrays. We have carried out Monte Carlo simulations, which 
agree with our experimental results. 
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I. Introduction 

Thermometers can be divided into two categories: pri- 
mary and secondary ones. The former do not have to be 
calibrated, but are generally slow and expensive, and are 
therefore mainly used for calibration purposes. Secondary 
thermometers are, on the other hand, usually cheap and 
fast, but they need to be calibrated against a known tem- 
perature scale, i.e. a primary thermometer or a scale which 
can be traced to one. The vast majority of thermometers 
used for routine measurements are secondary thermome- 
ters. 

Recently, a new primary thermometry method, the 
Coulomb blockade thermometry (CBT), was demonst- 
rated by Pekola et a/.|l]] It is based on the charging ef- 
fects in one-dimensional (ID) arrays of very small tunnel 
junctions. When the charging energy of the junctions, 
Eq — e 2 /2C, is greater than or of the order of the ther- 
mal energy, fc^T, the finite charge of the electron gives 
rise to a decreased differential electrical conductivity at 
low voltages (C is the tunnel junction capacitance and T 
is the temperature). This is called the Coulomb blockade 
of tunnelling ||. The Coulomb blockade shows up as 
a bell-shaped dip (figure |l|) when the differential conduc- 
tance, is plotted against the bias voltage, V. The 
depth of this dip is inversely proportional and the half- 
width is directly proportional to the temperature. Inter- 
estingly, the half-width depends only on the temperature, 
T, and natural constants, which makes it useful for pri- 
mary thermometry. This method has the advantage over 
other primary methods used at cryogenic temperatures in 
that it is both faster and cheaper, furthermore the physi- 
cal size of the thermometer can be made less than a mm. 
The depth of the conductance dip depends on the tun- 
nel junction capacitace, C, as well. Pekola et al. have 
demonstrated that this method can give an accuracy in 
temperature of at least 0.5% and they have used it suc- 
cessfully in the range from 20 mK to 30 K, more than three 



orders of magnitude, which is a large range for a primary 
thermometer. 

In this paper we show experimental data for a CBT con- 
sisting of a two-dimensional (2D) array of tunnel junctions, 
rather than the ID-arrays demonstrated by the group of 
Pekola et al. jl[ |[ [l2| One practical problem with ID- 
arrays of tunnel junctions is that they are dependent on 
every junction to function properly, like a chain is only 
as strong as its weakest link. So if one junction is dam- 
aged, ages rapidly or simply is bad from the beginning, 
the array might be unusable or at least its performance 
will be degraded. One way to reduce this degradation due 
to bad junctions is to connect several arrays in parallel, 
thus reducing the error to the order of 1/M, where M is 
the number of arrays. Even better is to make a 2D-array, 
where the metal islands are connected with tunnel junc- 
tions both in parallel and in series. The current would then 
simply pass around any broken junctions. This reduces 
the error due to a broken junction to the order 1/MN, 
where M and N are the number of junctions in parallel 
and in series respectively. (The resistance of an array of 
resistors increases with a factor of approximately 2/MN 
if one resistor is removed.) 

We have measured current-voltage characteristics on a 
2D array of tunnel junctions at temperatures from 1.5 K 
to 4.2 K in a pumped 4 He cryostat. Using the formulas de- 
rived for ID-arrays we have calculated the temperatures 
from these IV-curves and compared them with tempera- 
tures calculated from the pressure of the He gas in the 
cryostat. The correspondence is very good. 

We have also made Monte-Carlo simulations of 1D- 
and 2D-arrays and these simulations show that the tunnel 
junction capacitance C which is used in the equations for 
ID-arrays should be replaced in the 2D case by an effective 
capacitance C c ff which is higher than the real capacitance. 
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FIG. 1. The differential conductance as a function of bias voltage 
shows a bell-shaped dip around zero voltage. These measurements 
were made on a 256 X 256 junction 2D-array with a resistance of 
17kf2 and an effective capacitance of 2fF. The calculated curves 
were calculated with higher order corrections [Q. 

II. Theory 

The differential conductance G of a ID-array of N junc- 
tions has been calculated theoretically Q to first order 
in the limit Eq <C fegT. For arrays where all junctions 
have the same resistance and capacitance and negligible 
capacitance to ground, the result is 
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where the parameters un and vn are defined as 
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Gt is the differential conductance at voltages far above 
the Coulomb blockade. The effective charging energy, 
Ec,eS is defined similarly to the usual charging energy, 
except the physical capacitance is replaced with the ef- 
fective capacitance (which in a ID-array with negligible 
capacitance to ground is C c fr = jtzjC). 

The half-width Vl of the resulting conductance dip de- 
pends only on N and T: 



eVi = 5A39Nk B T. 



(4) 



The depth ^ = ^u N of the conductance dip also de- 
pends on the effective capacitance C e s of the junctions. 

At low temperatures (k B T w Eq) higher order correc- 
tions should be include d |4[ . This gives a correction to the 
halfwidth in equation (J4j), which is independent of tem- 
perature: 



eVi = 5A39Nk B T 1 + 0.3921 
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FIG. 2. Comparison between the temperature calculated from the 
half- width of the Coulomb blockade in the IV-curve and the temper- 
ature calculated from the 4 He equilibrium pressure. 



III. Sample parameters 

The measured array consisted of 256 x 256 tunnel junctions 
in a square pattern with a side of 200 /zm and the area of 
the overlap junctions was approximately 150 nm squared. 
The tunnel junctions were made of aluminium with a tun- 
nel barrier consisting of aluminium oxide, AlO^, fabricated 
with standard shadow evaporation technique J?], ||. The 
junctions had a tunnelling resistance of 17 kQ and a capaci- 
tance of 1 . 1 fF which means Eq /k B was . 85 K . The capac- 
itance was calculated from the offset voltage, V b =9.0 mV, 
using the offset analysis method described in ref ||. We 
have furthermore assumed that the offset voltage of a 2D- 
array of tunnel junctions is V g = eN/AC || . 



IV. Experimental setup 

The measurements were carried out in a pumped 4 He cryo- 
stat which has a temperature range from 4.2 K down to 
approximately 1.2 K. A custom built biasing and ampli- 
fication box fed the source signal to the sample and am- 
plified the current and voltage signals, which were mea- 
sured with two Keithley 2000 multimeters. A computer 
controlled Keithley 213 voltage source was used for the 
voltage sweep. The 4 He bath pressure was measured by 
a Wallace & Tiernan FA129 analog pressure meter which 
was connected directly to the pumped helium space. 



V. Results 

We measured IV-curves at 50 different temperature points 
and differentiated them numerically. Using equation |^ we 
calculated the temperature and compared it with the tem- 
perature calculated from the 4 He equilibrium pressure in 
the He bath, using the definition of ITS-90 ||. As can 
be seen in figure the two temperatures are almost the 
same, but not quite. The temperature from the IV-curve 
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FIG. 3. The relative difference between the temperature calculated 
from the IV-curve and from the 4 He pressure. There is a systematic 
difference of around 2%. As a comparison the uncertainty from the 
pressure reading is included, but it can't account for all the difference 
in temperature. 



is consistently slightly higher than the temperature from 
the 4 He pressure. 

The difference is clearer if we plot the relative difference, 
i.e. (Tjy — T pres )/T pres (figure ||). We see that there 
is a positive difference of 1-3% in all the measurements. 
The uncertainty in the pressure readings was estimated to 
about 1 Torr and the resulting temperature uncertainty 
is also shown in figure |[ While it may account for some 
of the spread, it can not explain the systematic difference 
of about 2%. One possible explaination could be that the 
IV-curves were measured quite rapidly one after the other, 
and perhaps the sample did not quite cool down at each 
temperature before the measurement was made. 

We have made Monte Carlo simulations of one- and two- 
dimensional tunnel junction arrays with a program pack- 
age called SIMON |HJ which is intended for simulations 
of circuits with very small tunnel junctions. Because this 
method turned out to be very time-consuming, we settled 
for quite small arrays, only three junctions in series and 
two in parallel. Both ID and 2D arrays behaved according 
to theory, except that in the 2D case the capacitance in 
equation ^ was not the actual capacitance of the tunnel 
junctions, but an effective capacitance C e g which was ap- 
proximately 1.4 times the actual value in the 3x2 junction 
array. This is not surprising, though, since two islands are 
not only connected via the junction between them, but 
also via the other islands in the array. In an infinite 2D 
array the capacitance between two neighbouring islands 
is exactly twice the capacitance of the tunnel junctions, 
and this is the reason that the offset voltage is a factor 
of two lower in a 2D-array compared to a ID-array [||. 
By calculating the effective capacitance from ^S- we get 
C e ff w 2.2 fF, which is twice the capacitance calculated 
from the offset voltage, as expected. 

An important property of the 2D-array is that it can 
be fabricated with lower resistance than a ID-array, even 
if it contains many junctions in series. This means that 
the measurement can be done faster and the measurement 
error is lower. When deducing the temperature from the 
half- width of the conductance vs voltage curve, the uncer- 
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FIG. 4. The simulated half-width voltage, scaled by the ideal value, 
plotted against A max /A m i n (dots) for a 5 X 5 junction array. j4 max 
and j4 m i n are the maximum and minimum junction areas respec- 
tively. The dashed lines are the corresponding results for ID-arrays 
by Hirvi et al. ]ll| , and the solid line is a quadratic fit to the simu- 
lated values, to compare with the quadratic dependence in the ID- 
array. It is clear that 2D-arrays are less affected by variations in the 
tunnel junction areas than ID-arrays, both at moderate (big figure) 
and large (inset) variations. Larger arrays give essentially the same 
result. 



tainty in T depends both on the uncertainty in G and in 
V. The measurement uncertainty can be written 
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where SI and SV are the current and voltage uncertain- 
ties respectively, given by the array resistance and tem- 
perature, and the amplifier noise. AI and AV are the 
intervals used for the numerical differentiation of the IV- 
curves or, if lock-in measurements are used, the modula- 
tion ampliture. Note that AI — GiAV where Gi is the 

2 2 

differential conductance at the half-width voltage. The 
last term in eq. ^| is the maximum differentiation error. 
Equation ^| has a minimum with respect to AV and we 
can chose AV (in our case approximately 5 % of Vi) to 
minimise the total error: 
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i?T is the tunnelling resistance per junction at the halv- 
width voltage. Note that the uncertainty decreases both 
with increasing N and M. 

Now we want to compare the measurement uncertainty 
of a ID- and a 2D-array. In this example we use a typical 
ID-array of 30 junctions (typical because this is the ap- 
proximate size used by Pekola et at), each junction with 
an effective capacitance of 2fF and a resistance of 20kf2. 



For the 2D-array we assume 256 x 256 junctions and the 
same junction parameters. The temperature is 4K. For 
the measurement we assume AD743 operational amplifiers 
at the input stage, because they have low voltage noise. 
The resulting uncertainty is 1.2- 10~ 3 for the ID-array and 
0.22 • 10~ 3 for the 2D-array, i.e. five times lower in the 2D 
case. One might argue that we would get the same effect 
if we fabricated 256 parallel ID-arrays with 256 junctions 
each, but then we might just as well make it 2D, for the 
sake of robustness, as discussed earlier. 

Another important issue is the tolerance to variations in 
the tunnel junction properties. Since no fabrication pro- 
cess is perfect, small variations can always be expected 
and in general, the smaller the structures, the bigger are 
the relative differences. Hirvi et al. have calculated the 
measurement error which can be expected in a ID-array 
at different inhomogenities [jllj . We have simulated a 5 x 5 
junction 2D-array with random distributions of tunnel 
junction areas and plotted the resulting half-width volt- 
age Vi against the parameter A max /A m - ln , where A max and 
j4 m i n are the maximum and minimum junction areas re- 
spectively (figure ||). Increasing the number of junctions 
give essentially the same result. We have assumed that 
RijCij — constant for all junctions, which follows from 
uniform tunnel barriers. It is clear from the figure that the 
2D-array is much less affected by variations in the junc- 
tion areas than ID-arrays. For example, at a 10% spread 
(A max /A m i„ = 1.10) the error in a ID-array is w0.15%, 
while in a square 2D-array it is «0.05%. It should be noted 
that even a 2 x iV-array is significantly better than a ID- 
array in this respect, at 10% spread the error is «0.08%. 

VI. Conclusions 

We have shown that a 2D-array of tunnel junctions can 
be used as a primary thermometer in the same manner as 
ID-arrays demonstrated earlier. 2D-arrays have certain 
advantages over ID-arrays for this application, namely 
the robustness to damage of individual tunnel junctions, 
higher tolerance to variations in tunnel junction proper- 
ties, lower signal noise and lower resistance which enables 
faster measurement. This should make them more suit- 
able for temperature measurements at temperatures above 
a few hundred mK where the electron heating effect is not 
a big problem. At lower temperatures ID-arrays are prob- 
ably better, since the topology of ID-arrays allow big cool- 
ing fins to be attached to the metal islands of the array, 
and this will improve the thermalisation of the electron 
system Q. 
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